GLOBAL ASYMPTOTICS OF STIELTJES-WIGERT POLYNOMIALS 
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Abstract. Asymptotic formulas are derived for the Stieltjes-Wigert polynomials S n (z;q) 
in the complex plane as the degree n grows to infinity. One formula holds in any disc 
centered at the origin, and the other holds outside any smaller disc centered at the origin; 
the two regions together cover the whole plane. In each region, the g-Airy function A q (z) is 
used as the approximant. For real x > 1/4, a limiting relation is also established between 
the g-Airy function A q (x) and the ordinary Airy function Ai(x) as q — > 1. 

1. INTRODUCTION 

We first fix some notations. Let k > be a fixed number and 
(1) g = exp{-(2fc 2 )- 1 }. 

Note that < q < 1. The g-shifted factorial is given by 

n-l 

(a,q) = l, (a,q) n = - aq 3 ), n = l,2,---. 

j=0 

With these notations, the Stieltjes-Wigert polynomials 

q 3 



(2) S n (x;q) = J2 * , (-xY, n = 0,l,2,---, 

{q,q)j{q; q) n -j 

are orthogonal with respect to the weight function 

(3) w(x) = krc'i exp{— k 2 log 2 x} 

for < x < oo; see (18.27.18)] and [8j (3.27.1)]. It should be mentioned that the Stieltjes- 
Wigert polynomials belong to the indetermined moment class and the weight function in ([3]) 
is not unique; see [3]. One important property of the Stieltjes-Wigert polynomials is the 
symmetry relation 

(4) S n {z-q) = {-zq n ) n S n (-^-,q 

which can be easily verified by changing the index j to n—j in the explicit expression given in 
(J2J). In some literatures, the variable x in is replaced by q^x\ see, for example, Szego [TT] . 
Chihara [2], and Wang and Wong [13]. The notation for the Stietjes-Wigert polynomials 
used in these literatures is 

(5) Pn(x) = S n (qzx;q). 

v(?; q)n 



The work of Y. T. Li is supported by the HKBU Strategic Development Fund and a start-up fund from 
HKBU (Fund Ref.: RC-NAC AD-LI YT). 

1 



2 



Y. T. LI AND R. WONG 



The Stietjes-Wigert polynomials appear in random walks and random matrix formulation 
of Chern-Simons theory on Seifert manifolds; see [UH]. 

The asymptotics of the Stieltjes-Wigert polynomials, as the degree tends to infinity, has 
been studied by several authors. In 1923, Wigert [T5] proved that the polynomials have the 
limiting behavior 

1/4 00 k 2 +k/2 

(6) hm (-ir q ^ 2 Pn ( X ) = - 7 L== j2(-v k V^T**' 



which can be put in terms of the g-Airy function 

CO fe2 

This function appeared in the third identity on p. 57 of Ramanujan's "Lost Notebook" [10J. 
(For this reason, it is also known as the Ramanujan function.) In terms of the g-Airy 
function, Wigert's result can be stated as 

(8) lim S n (x;q) = \ A q (x). 

It is known that all zeros of S n (x;q) lie in the interval (0, 4g _2n ); see [Ti] . Hence, we 
introduce a new scale 

(9) z := q- nt u 

with u G C \ {0} and t Gt. The values of t = and t = 2 can be regarded as the turning 
points of S n (q~ nt u; q). Taking into account the symmetry relation in (J4J), one may restrict 
oneself to the case t > 1; see [T31 (1-4)]. (However, in the present paper, we will not make 
this restriction.) The case t = 2 has been studied by Ismail [6], and he proved 

(10) lim q n \-u) n S n {uq- nt -q) = - — -r — A q (-), t = 2, 

uniformly on compact subsects of C \ {0}; see Theorem 2.5]. This result can in fact 
be derived directly from Wigert's result in ([8]) via the symmetry relation mentioned in (j4]). 
In [7], Ismail and Zhang extended the validity of this result to t > 2. For 1 < t < 2, Ismail 
and Zhang [7] gave asymptotic formulas for these polynomials in terms of the theta-type 
function 



but in a very complicated manner. The result in [7] was then simplified by Wang and 
Wong 1121 . For instance, when 1 < t < 2, Wang and Wong proved that 



( _ r ,,\n-m n n 2 (l-t)~m[n(2-t)-m] ( / n 2m-n(2-t)\ ~\ 

(12) S ^"' 5) = - — { e « (^r~) + ^ } • 

where / = |(2 — t), m = \nl\ and 5 > is any small number; see [121 Corollary 2]. Note 
that all these results are not valid in a neighborhood containing t = 2, one of the turning 
points. To resolve this issue, a uniform asymptotic formula was given by Wang and Wong in 
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a second paper [13J. For z :- 
they showed that 



uq 



-nt 



with t > 2(1 — 8), S being any small positive constant, 



(13) 



S n (z;q) 



-z) n q n 
(<?; q)n 



[A q , n (q~ 2n /z) + r n (z)}, 



where r n (z) is the remainder and A q ^ n (z) is the g-Airy polynomial obtained by truncating 
the infinite series in ([?]) at k — n, i.e., 



k=0 



i-z) k . 



In this paper, we shall show that the g-Airy polynomial A qtU (z) in equation (|T3|) can be 
replaced by the g-Airy function A q (z). Moreover, we shall show that the resulting formula 
is global. More precisely, we have the following result. 



Theorem 1. Let z := uq with —oo < t < 2, u G C and \u\ < R, where R > is any 
fixed positive number. We have 

(14) S n (z;q) 

where the remainder satisfies 



[A q (z) +r n (z)] 



(15) 



Vn(z)\ < 



q 



n(l-cr) 



+ 



1 — q 1 — q 



[no-J 



A q {-\z\) 



with a = max{i, ^ + |}. 



-2' 2 1 4- 

Let z := uq~ nt with < t < oo, hgC and \u\ > 1/R, where R > is any fixed positive 
number. We have 

,2 



(16) S n (z;q) 
where the remainder satisfies 



-z) n q r 



(q; q)i 



[Mq- 2n /z) + 



r n (z 



(17) 



rJz) < 



q 



n(S-l) 



l-q 1 - q \2 



[n(2-5)\ 



A q (-q- 2n /\z\) 



with 5 = min{|, 1 + |}. 

Note that the quantities inside the square brackets in (1151) and ( II 7p are exponentially small. 
Furthermore, the sizes of the g-Airy functions in these two equations for large values of their 
arguments are about the same as the leading terms in their corresponding approximation 
formulas in (IHD and flT£); cf. (£2} below. 

The investigations mentioned above all started with the explicit expression of S n (x;q) 
given in 02]). In [13], Wang and Wong used a different method, namely, the Riemann-Hilbert 
approach, to get a uniform asymptotic expansion of the Stietjes-Wigert polynomials p n (x) in 
terms of Airy functions. But, the main result in [T3] needs a correction, and the correction 
is that the parameter k in ([[]) and (EJ) should depend on n and tend to infinity as n tends to 
infinity. In other words, the result in [13] holds with a varying weight and the number q in ([T]) 
is required to approach 1. In fact, more precisely, if k ~ n a as n — > oo and < a < |, then 
all formulas in [T3] remain valid; if a > \, then some equations need be amended, but the 
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main result still holds; see also Baik and Suidan [Tj. Note that the results in Theorem [T] hold 
even when q tends to 1. Comparing the results in Theorem [1] and in |14j . we will establish 
the following limiting relation between the g-Airy function and the ordinary Airy function: 



Theorem 2. Let £(x) denote the function defined by 

ris) 



2 l /-l°g(4x) 

- [£( x )] 3 ^ 2 = "i — 7 — rr / arctan \/ e s — Ids. 

3 log(l/g) Jo 

Then, for fixed x > 1/4, the q-Airy function has the following asymptotic approximation 



(19) 

as q —7- 1. 



A q {y/qx) ~ 2A/7rexp 



3 log X — 7T 2 



121og(l/g) J \Ax-l 



Ai(-£(s)) 



This result in fact holds for any complex x, if we replace £(x) in (fl~9l) by defined by 

1 



lo s(4z) , \ofr 2 (4r) 

log (l + VT^) ds - pp 



3 L J log g 

In view of this result, A q is indeed a g-analogue of the Airy function. 

2. Global asymptotics of Stieltjes-Wigert polynomials 

Proof of Theorem 1. 

If u = 0, then z = uq~ nt = and 

(g;g) n S' n (0;g)=^(0) = l 

by the series representations in (TS]) and (j7|). Hence, the result in (|14p - (|15p follows immediately, 
and we just need consider the case when u^O. Noting that 

-nt U - n (t-=*W) 

z = uq =- — -q v «iog9/ ) 
\u\ 

we may assume \u\ — 1 without loss of generality 
For notational convenience, we put 

r n{z) =(q; q) n S n (z; q) - A q {z) 



(20) 



where 



E 

j=0 



(9;?)r 



.(?5?)n-i 

=:/i + J 2 + J 3 , 



- 1 



9 



(?; ?)j 



-2 



-z 



[na\ 
3=0 



E 

j=|_ncrj+l 



(g; gQn _ 
{q;q)n 



(?; 9)n- 



- 1 
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and < cr < 1 is a constant to be specified later. In view of the inequality 1 — ab < 
(1 — a) + (1 — b) for any a, b G (0, 1), we have 



(21) 1 - = 1 - (q^; q) 3 < £ g -*H < I . 

Thus, we have 



n-j 



(22) Jx < £ g g k < ^ A q {-\z\). 

Let 

f 1 t 1 1 f n , ni . 1 

(23) a = max < -, — I — > , / = max < — , — + 1 > . 

It is readily seen that | < a < 1 and that y < / < n for sufficiently large n. For j > |_?7.crJ 
and n sufficiently large, there exists 7 > such that 

(j 2 - ntj) - (/ 2 - ntl) = (j + nt) > (j - I) 2 > 7 2 J 2 . 

Then, it follows that 



OO ,-2 

q 3 



\h+h\<-±- £ ■ , 



z\ 



j=\n<r\+l 



ql 2 -ntl 00 q(j 2 -ntj)~(l 2 -ntl) 

< — v <? 7V . 



Since < | for j > [ricrj and sufficiently large n, we have 

lr„ 4. r„l <. 

(1 - g)(g; g) 



J 2 -ntl 00 

j 2 + j 3 i< ? y ' x 



^24) "' ■ °° i=L no "J+ 1 

I 2 -ntl /^\ L«°"J 

< — 

(1 - g)(g; g) c 
Similar to the inequality in (I2ip . we have 

i-^% = i-(g m ;g)oo<^<^<- 

(q;q)i l-q 1-q 2 

for sufficiently large n, which leads to 

(25) 7 ^<^ T <2. 
Moreover, we have 

l 2 -ntl 00 n j 2 -ntj 00 j 2 

(26) ±— — < V -— — < V \q- nt u\ j = 

{q;q)i (q;q) 3 
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A combination of the inequalities in fl24")) . ( |25|) and fl26|) gives 

2 /\ \ L n<J J 



which, together with p)} and ([22]), yields f lT5]) . 

The result in f JT6|) - ffT7|) can be proved in a similar manner. One can also obtain this directly 
from f fl4|) . ffT5|) and the symmetry relation of ^(z; g) mentioned in (jl]). 

3. Limiting behavior of A q (z) as c? ->■ 1 

It is known that A 9 (z) has infinitely many positive zeros and satisfies the three-term 
recurrence relation [5] 

(27) A q {z) - A q (qz) + qzA q (q 2 z) = 0. 
Moreover, Zhang [16] has shown that 

(28) lim AMl -q)z) = e~ z 

for any fixed z G C. In [T3J, Propsition 1], Wang and Wong proved that 

(29) A q (z) = K ~ Z) . g [e,(-g 2m ^) + (q m ^) } 

as z oo, where m := j and 5 > is any small number; see also Zhang [Ml Theorem 

2.1]. 

In this section, we shall establish the limiting relation of A q (z), as q — > 1, stated in 
Theorem [2J Let us first review some of the results given in Wang and Wong [14J . Let 

(30) k:=k{n)=n^. 
Then 

(31) q = exp{-(2fc 2 )- 1 } = exp{-(2 v ^)- 1 } = l--L= + o(^j 

as n — > oo. In Sec. 1 we have commented that the main result in Wang and Wong [13] holds 
if the parameter A; in (TTJ) depends on n and satisfies a growth condition such as the one given 
in ( 1301) . The MRS numbers a n and (3 n have been calculated in [H], and are explicitly given 
in equations (3.19) and (3.20). We note that 

(32) ^ ~ ^q~ (2n+1 \ asw^oo; 

cf. [IH eq. (2.4)]. Furthermore, the change of variable t —¥ y defined by 

t 



(33) y = a/ a n (3 r , 



cxp 



2 log(/3 n /a n ) 



takes the interval — 1 < t < 1 onto the interval a n < y < (3 n , where y = l/(xq 2n+1 ); see [HI 
(2.14)]. Let 7r n (x) denote the monic Stieltjes-Wigert polynomial 

/ \ Vn{x) 
7T n (x) 



In 
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where p n (x) is the polynomial given in and 

7n = g nW/ 7v / Mk 

that is 

(34) n n (x) = (-l)«g-" 2 -«/ 2 (g; q^S^q^x; q). 
The major results in [H] is the asymptotic formula 

(35) n n (y) = {ivUi A(y, n) + O (n~$) 



where iV = n + 4 



2' 



f36) A(v n) - Z ^ 

(36) AM ~ 1(V -a n )(y -AOF 4 ' 



and 



-[— ?? n (t)l 3 ^ 2 = --r-; — t- . . / arctan ^ ^ -rfr 

3 L /nv yj N\og{\/q) J t e aT + 1 

« : = o log — ; 

2 a„. 



see [TU (2.11), (2.12) and (6.7)]. Recall that y in fl35|) and fl36l) is a function of £, given in 
033]). Using the formula [H (6.8)] 

|-"" (t)13/2 = m4m f " arctan + ° ( al "0 ' 

where — 1 + S < t < 1, we have 

(37) ivi^ ~ 

as n — >■ oo, where is the function defined in (ITS]) . The last equation gives 

(38) Ai(N 2/3 r] n ) ~ Ai(-f(z)), n ->■ oo. 
Next, we show that 

(39) ivU(y,n)~2v^^^Ly , n -> oo. 
To this end, we note from (l36l) that 



(40) 



' ' xg 2 ™ +1 ' J \[l/xq 2n+1 -a n ][l/xq 2n+l - (3 n \ 



\ N 2 ISn{Wn-an) 2 X 2 q 2{2n+1) ^'' 

\{l-a n xq 2n + l ){l-^ n xq 2 ^) 
By using equations (3.17) and (3.18) in [TJ], we have 
(41) a n f3 n = q~^ +1 \ 

Since (3 n is large, the last two equations, together with (|37"j) . give 



X q2n+i i J \ 1 — x/a 
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thus proving (1391) . Here, we have also made use of the fact that a n ~ \ as n — > oo. 
Finally, we evaluate the asymptotics of 

/.ox \Re^ ln ,- fl, 1,,, 2 1, k , 1 

(43) == = Vvrexp^ -ln + ^k log y--\og—=-n\ogy). 

Recall y = l/xq 2n+l . Using the formula [H (3.29)] 

MAT-1) k 2 7r 2 , fc , Ar 

we obtain 

(44) ->^= ~ v^Fexp { N{N ~ 2 l) - k -f+ l -k 2 \o g 2 {xq 2 ^)+n\o g {xq 2 ^)\. 

Since log(l/g) = l/2k 2 and iV = n + |, one can show that the right-hand side of (144j) is 
equal to 

/ 1 3lQg 2 X-7T 2 | /¥ f3l0g 2 X-7T 2 

(45) V^exp ^ — o io & x + "TTTi — / \ f = U ~ ex P 



2 & 121og(l/g) j I 121og(l/g) 

Hence, we obtain 

y/^Fe^" /tt / 3 log 2 g - 7T 2 | 

(4dj . ~ \ / — exp < — ; — — — > , n — > oo. 

1 ' ^V 7 ^) V * I 121og(l/g) J 



A combination of ( 1351) . (1391) and (1461) yields 

(47) ( V»«)".„(l/ I5 -Vyfe X p{^f^}2^(^) i Ai(- fW ). 
Therefore, we have 

(xq 2n+1 ) n 7c n (l/xq 2n+1 ) = (q;q)n(-q 1/2 xq n ) n S n (q- 2n /q 1/2 x;q) = (q; q) n S n (q*x;q). 

Here, we have made use of the symmetry relation of S n (z;q) given in (jj]). Note that the 
result in Theorem [1] in fact holds with q satisfying (I3TI) . Thus, Theorem [1] gives 

(48) (xq 2n+1 ) n 7r n (l/xq 2n+1 ) = A q (^x) + r n (y/qx). 
Coupling ( )47l) and ( |48l) . we obtain the desired formula 

a / / — \ n r- f 31og 2 x-vr 2> l / £{x) \* ^ v, 

as g — t- 1~. 

4. Numerical Verification 

In the following tables, we have used a Maple-aided program to verify the asymptotic 
formulas in Theorem [1] and the limiting relation in Theorem [2j The values are represented 
in scientific notations; for example, 

-2.325e-3 = -2.325 x KT 3 = -0.00235. 

In Table 1, 'True' stands for the true value of (q; q) n S n (uq~ nt ; q) by summing the series in 
d2J); 'Approx.' stands for the approximate value obtained from (fT4"|) ; 'Error' is the relative 
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error of the approximate value. The degree of the polynomial is n = 50 and q = 0.5 in Table 
1, u takes the values of 1, —1 and 1 + i, and t takes the values of 0, 0.5, 1.0, 1.2 and 1.6. Here, 
we would like to mention that the approximate values are very close to the true values. For 
instance, in the case when u = t = 1, the true value is -5.83981318477869- • • x 10 187 , whereas 
the approximate value is -5.83981318477868- • • x 10 187 . If we take only a few digits as what 
we have done in Table 1, the two values all appear to be nearly the same. In Table 2, 'True' 
stands for the true value of A q (^/qx) by summing the series in ((7j); 'Approx.' stands for the 
approximate value obtained from the quantity on the right-hand side of f JT9|) ; 'Error' is the 
relative error of the approximate value. We examine the cases when x = 0.5, 1.0,4.0, 10,20 
and q = 0.9, 0.92, 0.94, 0.96, 0.98, 0.99. 

Table 1 . Numerical Verification of Theorem 1 



\ t 

u \ 







0.5 


0.8 


1.0 


1.2 


1.6 




True 


0.16076 - 


9.3534e42 


1.0831el20 


-5.8398el87 


3.5453e270 


1.8649e481 


1 


Approx. 


0.16076 - 


9.3534e42 


1.0831el20 


-5.8398el87 


3.5453e270 


1.8649e481 




Error 


2.99e-15 


2.98e-8 


1.78e-15 


1.18e-15 


6.05e-13 


6.36e-7 




True 


2.17267 i 


3.0063e47 


1.9036el21 


1.0264el89 


6.2313e271 


3.2740e482 


-1 


Approx. 


2.17267 I 


3.0063e47 


1.9306el21 


1.0264el89 


6.2312e271 


3.2778e482 




Error 


6.31e-16 


6.12e-12 


l.lle-9 


3.54e-8 


1.13e-6 


1.16e-3 




\ t 

u \ 









0.5 


1.0 




1.6 










True 


0.0117-0.6786i (1.92+8.38i)e48 


(-8.18-1.87i)el91 (4.107-2.571i)e487 






1+i 


Approx. 


0.0117-0.6786i (1.92 


+8.38i)e48 


(-8.18-1.87i)el91 (4.106-2.578i)e487 










Error 


8.83e-17 3.16e-12 


1.39e-8 


4.55e-4 










Table 2. Numerical Verificat 


ion of Theorem 2 






\q 

X \ 




0.9 


0.92 


0.94 


0.96 


0.98 


0.99 










True 


-2.325e-3 


-3.826e-4 


1.120e-5 


-2.966e-8 


5.080e-16 


4.9298e-32 






0.5 


Approx. 


-2.320e-3 


-3.819e-4 


1.118e-5 


-2.964e-8 


5.078e-16 


4.9303e-32 










Error 


0.0022 


0.0018 


0.0012 


0.00080 


0.00038 


0.0000995 










True 


-5.171e-4 


2.978e-5 


-2.556e-6 


1.326e-9 


2.178e-18 


4.1417e-36 






1.0 


Approx. 


-5.159e-4 


2.973e-5 


-2.553e-6 


1.325e-9 


2.177e-18 


4.1408e-36 










Error 


0.0022 


0.0018 


0.0013 


0.00087 


0.00043 


0.00021 










True 


0.034973 


0.01680 


4.4202e-4 


-1.0084e-4 


4.4912e-8 


-5.6869e-16 






4.0 


Approx. 


0.034891 


0.01677 


4.1898e-4 


-1.0073e-4 


4.4893e-8 


-5.6853e-16 










Error 


0.00235 


0.0018 


0.0028 


0.00107 


0.00043 


0.00028 










True 


38.6522 


-247.876 


2715.83 


43744.8 


3.3978el0 


2.1941e21 
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Approx. 


38.5316 


-247.372 


2712.29 


43745.2 


3.3961el0 


2.1944e21 










Error 


0.0031 


0.0020 


0.0013 


0.00022 


0.00051 


0.00014 










True 


-2.0951e5 


-3.5927e6 


-5.9716e9 


7.3472el4 


-6.1129e29 


1.5900e61 






20 




Approx. 


-2.0884e5 


-3.5801e6 


-5.9645e9 


7.3418el4 


-6.1124e29 


1.5897e61 










Error 


0.00320 


0.00349 


0.00119 


0.00073 


0.00007 


0.00023 
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